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Abstract. Radio emission models of colliding wind binaries (CWBs) have been discussed by Dougherty et al. 
(2003). We extend these models by considering the temporal and spatial evolution of the energy distribution of 
^ . relativistic electrons as they advect downstream from their shock acceleration site. The energy spectrum evolves 

^Nf^ ' significantly due to the strength of inverse-Compton (IC) cooling in these systems, and a full numerical evaluation 

QQ , of the synchrotron emission and absorption coefficients is made. We have demonstrated that the geometry of the 

CNj ■ WCR and the streamlines of the flow within it lead to a spatially dependent break frequency in the synchrotron 

' emission. We therefore do not observe a single, sharp break in the synchrotron spectrum integrated over the WCR, 

but rather a steepening of the synchrotron spectrum towards higher frequencies. We also observe that emission 
from the wind-collision region (WCR) may appear brightest near the shocks, since the impact of IC cooling on 
the non-thermal electron distribution is greatest near the contact discontinuity (CD), and demonstrate that the 
. impact of IC cooling on the observed radio emission increases significantly with decreasing binary separation. We 

CL(' study how the synchrotron emission changes in response to departures from equipartition, and investigate how 

JL I the thermal flux from the WCR varies with binary separation. Since the emission from the WCR is optically thin, 

^ ■ we see a substantial fraction of this emission at certain viewing angles, and we show that the thermal emission 

from a CWB can mimic a thermal plus non-thermal composite spectrum if the thermal emission from the WCR 
, becomes comparable to that from the unshocked winds. We demonstrate that the observed synchrotron emission 

depends upon the viewing angle and the wind-momentum ratio, and find that the observed synchrotron emission 
decreases as the viewing angle moves through the WCR from the WR shock to the O shock. We obtain a number 
of insights relevant to models of closer systems such as WR 140. Finally, we apply our new models to the very wide 
5^ ' system WR 147. The acceleration of non-thermal electrons appears to be very efficient in our models of WR 147, 

and we suggest that the shock structure may be modified by feedback from the accelerated particles. 

Key words. stars:binaries:general - stars;early-type - stars:individual:WR 147 - stars:Wolf-Rayet - stars: winds, 
outfiows - radio continuum:stars 



1. Introduction 

Radio observations of carly-typc stars reveal that they 
can be sources of both thermal and non-thermal radia- 
tion. The thermal emission is readily explained as free- 
free emission within the stellar wind (Wright & Barlow 
1975), while the non-thermal emission is attributed to rel- 
ativistic electrons moving in a magnetic field. The consen- 
sus is that the relativistic electrons in the winds of mas- 
sive stars are created through diffusive shock acceleration 
(DSA), though other mechanisms have been investigated 
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(e.g., Jardine et al. 1996). Observations of the WR+OB 
binary systems WR140 (Dougherty et al. 2005), WR146 
(Dougherty et al. 1996, 2000) and WR147 (Moran et al. 
1989; ChurchweU et al. 1992; Wilhams et al. 1997; Niemela 
et al. 1998) clearly show the non-thermal emission in 
these binaries arises from a wind-collision region (WCR) 
bounded by strong shocks formed where the two stellar 
winds collide. There is strong evidence that all WR stars 
that exhibit non-thermal emission are binary (Dougherty 
& Wilhams 2000). 

For synchrotron emission to be observed, the relativis- 
tic electrons must be within the optically thin region of the 
stellar wind(s), which in the case of massive stars implies 
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at more than a few hundred stellar radii from the underly- 
ing star. These large radii are problematic for single star 
models of synchrotron emission, where the declining ve- 
locity jump and compression ratio with increasing radius 
of wind-embedded shocks produces a rapidly decreasing 
synchrotron emissivity (Van Loo 2005). 

CWB systems are important objects for investigating 
the underlying physics of DSA because they provide ac- 
cess to higher mass, radiation and magnetic field energy 
densities than is possible through studies of supernova 
remnants. However, until the work of Dougherty et al. 
(2003) (hereafter Paper I), CWB models of radio spectra 
were typically based on a highly simplified model consist- 
ing of a thermal source plus a point-like source of non- 
thermal emission attenuated by free-free absorption (e.g.. 
Chapman et al. 1999; Monnier et al. 2002). Use of these 
models often resulted in difficulties in fitting observational 
data (see, e.g., WiUiams et al. 1990; White & Becker 1995). 
As a first step toward the construction of more realistic 
models with better predictive qualities, a hydrodynami- 
cal model of the WCR was used in Paper I to obtain a 
more accurate representation of the spatial distribution 
of the free-free and non-thermal emission from CWB sys- 
tems. Several important scaling laws were determined and, 
in spite of other simplifying assumptions, the results of 
that work showed great promise in comparisons with ob- 
served continuum spectra and images. We refer the reader 
to Paper I for essential background to the work presented 
here. 

Inverse-Compton (hereafter IC) cooling is an impor- 
tant energy-loss mechanism for relativistic electrons in 
CWB systems, since the energy distribution will evolve as 
the electrons advect away from the shocks. The inclusion 
of this process and the calculation of emission and absorp- 
tion coefficients from an arbitrary energy distribution of 
non-thermal electrons are the main enhancements to the 
model and the work presented in this paper. Following 
examination of the effect of these mechanisms, we apply 
these new models to WR 147, and close by describing fu- 
ture directions for this work. 

2. Modelling the radio emission from CWBs 

In Paper I. the radio emission and absorption from a 
CWB system was calculated from a 2D axis-symmetric 
hydrodynamical model of the stellar winds and the col- 
lision zone. The temperature and density values on the 
hydrodynamic grid were used to calculate the free-free 
emission and absorption coefficients from each grid cell. 
The synchrotron emission and self-absorption from each 
cell within the WCR was calculated assuming that the 
distribution of relativistic electrons in each cell could be 
specified by a power-law, i.e. Uei'y) (x 7~^, where 7 is the 
Lorentz factor^. The synchrotron emissivity, P^, also fol- 
lowed a power-law (Pi, cx i/"^^"-'^)/^), which was assumed 

^ For a strong shock with a compression ratio of 4, test par- 
ticle theory predicts that p — 2 (see references in Paper I). 



to be valid over all frequency, v. The fiux and intensity 
distribution at a specified frequency were then obtained 
from a radiative transfer calculation. 

The hydrodynamical simulations of Paper I do not pro- 
vide direct information on the magnetic field or relativistic 
particle distribution, so the magnetic energy density Ub, 
and the relativistic electron energy density U^ei, were set 
proportional to the thermal particle internal energy den- 
sity, Uth- 

Ub = B^Stt = CeC/th, (1) 
and 

Urel = J no(7)7TOoC^(i7 = CrelC^th- (2) 

C^th = r~L' '^itl^ P the gas pressure and F the adiabatic 
index, assumed to be 5/3, as for an ideal gas. Cb and 
Crci are constants which determine the amount of syn- 
chrotron emission and absorption, and the Razin turnover 
frequency in our models. Since we did not attempt to 
model the acceleration of the relativistic particles from ba- 
sic principles, the values of Cb and Crci are not determined 
in a self-consistent fashion. Instead, when modelling spe- 
cific systems, their values arc chosen to best match the 
observed radio emission. Our previous models in Paper I 
showed that if Cb = Crci, then values of around 1% are re- 
quired (sec also Mioduszewski et al. 2001). Unless other- 
wise noted, Cb = Crci ~ Ci though we relax this restriction 
in Sec. 3.6 and in some models in Sec. 4. 

2.1. The hydrodynamical code 

In contrast to Paper I, a hydrodynamical code with adap- 
tive mesh refinement is used here (Falle & Komissarov 
1996; Falle et al. 1998). Such codes concentrate more cells 
where the flow is interesting, for example near shocks, and 
fewer cells where it is uniform, resulting in significantly 
shorter computational times. We use a small amount of 
artificial viscosity in the calculations to prevent the de- 
velopment of the carbuncle instability at the apex of the 
WCR. This numerical artifact develops when shocks are 
locally aligned with the grid (Walder 1994; Levcque 1998). 

2.2. The non-thermal electron energy spectrum 

In this work, it is assumed that the acceleration of rel- 
ativistic electrons at the shocks bounding the WCR pro- 
duces an energy distribution with a power-law spectrum. 
When particle acceleration is relatively inefficient, this is 
only strictly true in the absence of cooling. However, the 
resulting deviation from a pure power-law is restricted to 
energies which are very near to the high energy cutoff 
(Chen & White 1991; Van Loo et al. 2005), and we ig- 
nore this effect in the current work. In contrast, when 
DSA is efficient, non-linear effects modify the pre-shock 
wind speeds over a length scale comparable to the diffu- 
sion length of the energetic particles, and slow them before 
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they undergo a shock transition. Thus higher energy par- 
ticles with larger diffusion length scatter in regions where 
the pre-shock wind speed is closer to the terminal speed of 
the stellar wind, while low energy particles do not diffuse 
very far from the subshock and scatter in regions where 
the wind speed is lower due to its journey through much of 
the shock precursor. This speed profile results in a concave 
upward curvature to the non-thermal particle spectrum, 
where p decreases with increasing energy (Ellison et al. 
2004, and references therein). Efficient particle accelera- 
tion can also flatten the non-thermal particle energy dis- 
tribution i.e. p < 2. Modelling this process is non-trivial 
and we make no attempt in the current work to calculate 
these effects self-consistently. However, we discuss it fur- 
ther and perform calculations with a non-thermal particle 
spectrum where p < 2 in Sec. 4. 

Relativistic electrons that scatter out of the shock with 
Lorcntz factor 7 < 10^ are frozen into the post-shock 
flow (White 1985). As they arc advected out of the sys- 
tem, the relativistic electrons lose energy through the ac- 
tion of several processes, including ionic cooling, magneto- 
bremsstrahlung cooling, and IC cooling. Since IC cool- 
ing strongly dominates magneto-bremsstrahlung cooling 
in this work (see Eq. 14 of Paper I), we do not consider 
the latter here. Cooling reduces the non-thermal electron 
energy density at a given 7 and alters the energy spec- 
trum of these electrons from a simple power-law. Thus, 
in this work we assume the non-thermal electron energy 
spectrum is a power-law at the shocks, but let it evolve 
as the electrons advcct downstream from the shocks. As 
noted earlier, the inclusion of these processes is the main 
enhancement to the model since Paper I, and in the fol- 
lowing sub-sections we describe in detail our method of 
implementation. 

2.2.1. Inverse Compton cooling 

The high UV luminosity of the stars in CWBs means that 
IC cooling of the relativistic electrons in the WCR can 
be very efficient, particularly for high 7 electrons. In the 
Thomson limit, the rate of energy loss of a relativistic 
electron exposed to an isotropic distribution of photons 
by IC scattering is (Rybicki & Lightman 1979) 



compt — T^OtC^I P J7ph, 



(3) 



where ut is the Thomson cross section, (3 = v/c = 
•\/l — 1/7^, U-pYi is the energy density of photons, and 
where elastic scattering is assumed and quantum effects 
on the cross section arc neglected. For a single star and 
relativistic particles with /3 w 1, we can rewrite Eq. 3 as 



(i7 



(7^7 



L 



= 8.61 X lO^^o^Vs"^ 



(4) 



where L is the stellar luminosity, r is the distance of the 



particles from the star, and / 



( erg s ^ cm ^ ) is a 



For a given electron, the reduction in 7 due to IC cooling 
can be calculated if the integral of ft along its streamline 
is known. In a CWB, / = + due to the radiation 
field of both stars. The value of ft in any given cell on 
the hydrodynamic grid can be obtained by the specifica- 
tion of a suitable scalar in the code. It is then a simple 
matter during post-processing to calculate the IC-coolcd 
relativistic energy spectrum at that location. 

The influence of IC cooling increases in magnitude as 
plasma flows from the shocks toward the contact discon- 
tinuity - this reflects the increased exposure to the high 
intensity stellar radiation field from the time when the 
relativistic electrons were accelerated in one of the shocks 
bounding the WCR. In addition, since the rate of cooling 
is proportional to 7^, IC cooling is most significant for 
electrons with the highest energies. We can therefore de- 
fine a characteristic Lorentz factor, 7c, where IC cooling 
becomes significant for flow near the apex of the WCR. 
The flow morphology and the distribution of electrons 
with a specific Lorentz factor is shown in Fig. 1. For low- 
energy electrons (7 << 7c), IC cooling is neghgible and 
such electrons may exist throughout the WCR, as shown 
in Fig. lb. In contrast, IC cooling is significant for elec- 
trons with 7 > 7c, which rapidly lose energy as they flow 
away from the shocks. In such circumstances, electrons 
with 7 ^ 7c are effectively conflned to narrow regions 
near the shocks, and absent from the central volume of 
the WCR, as illustrated in Fig. Ic. 

2.2.2. Ionic (coulomb) cooling 

Relativistic electrons may also lose energy through 
Coulomb collisions with thermal ions. The rate of energy 
loss through this process for electrons of energy 7 is (cf. 
Chen & White 1991) 



(i7 



coul 



47re^ni In A 



1 



3.0 X 10 



— 1/2 

where In A — 29.7 + In n-^ Tq is the Coulomb loga- 
rithm (e.g., Cowie & McKee 1977)^, Tg is the temperature 
in units of 10^ K, and rii is the number density of ions 
( cm"'^ ). As for IC cooling, we can estimate the degree of 
ionic cooling, namely by constructing a scalar to track the 
value of Tiit In A in each hydrodynamic cell of hot shocked 
gas. Ionic coooling is greatest for non-thermal electrons 
with relatively low energies, and while it is included in our 
model we note that there is no observational signature of 
this process on the synchrotron spectrum above 10 MHz 
(see Fig 3). However, we include it for completeness. 



_X4 n-i In A _^ 



(5) 



2.2.3. Maximum energy of relativistic electrons 

In the standard DSA scheme, particles arc assumed to 
scatter from MHD waves in the background plasma, either 

^ While we use this expression for the Coulomb logarithm, 
Li (2001) note that it should be reduced to about half of its 



measure of the radiative flux on an electron at distance r. conventional value. 




Fig. 1. Schematic diagrams of the fiow within the WCR and the distribution of relativistic electrons, a) The flow within the 
winds and the WCR, which is located nearer the star with the weaker wind (star 2 in this example). The dotted line depicts 
the contact discontinuity (CD) between the two winds. Either side of the CD is a shock which separates the high Mach number 
wind from high pressure post-shock gas. Example streamlines are shown for material leaving the stars along specific direction 
vectors, b) The distribution of relativistic electrons with Lorentz factor 7 << 7c, assuming that such electrons are first-order 
Fermi accelerated at the shocks, c) As b) but for 7 '-^ 7c. IC cooling effectively excludes electrons with Lorentz factor 7c from a 
region near the CD. 



pre-existing or generated by the counter-streaming ions 
themselves (Bell 1978). There is strong evidence that the 
efficiency of particle acceleration in the DSA scheme is 
dependent on the orientation of the upstream magnetic 
field to the shock normal (see the review by Blandford & 
Eichler 1987). When the shock is quasi-perpendicular the 
competing "shock drift" mechanism may be more efficient 
(Jokipii 1987), though the acceleration rates become equal 
in shocks with a high level of magnetic turbulence since 
these have no well-defined obliquity'^. 

Since the physics of electron injection is not yet fully 
understood, we ignore these complications in the present 
work, and assume that both the normalization of the non- 
thermal particle spectrum, and the maximum energy that 
the accelerated particles can reach, are independent of spa- 
tial position. To estimate the value of 7max, we balance the 
rate of energy gain through Ist-order Fermi acceleration, 
by the rate of energy loss due to IC cooling (the dominant 
cooling mechanism for relativistic electrons with large 7). 
The Ist-order Fermi power is (Chen & White 1991; Van 
Loo et al. 2005)^ 



It 



4 Xc - 1 ti7/3^ 2ii 
3' 



A 



(6) 



where Xc is the ratio of upstream to downstream flow 
speeds, u is the shock velocity, c is the speed of fight, 
and A is the mean-free path of the accelerating charge. 
This latter quantity is somewhat uncertain, but is typi- 
cally estimated as several gyro-radii. We adopt A « 3rg, 
where the gyro-radius rg = pc/qB (White 1985)^, and p 



^ The obliquity of the shocks expected in CWB systems is 
described by Eichler & Usov (1993). 

* Chen & White (1991) erroneously assume that shock- 
accelerated electrons return to the shock after one scattering, 
whereas electrons are actually scattered c/2u times before re- 
crossing the shock (Van Loo et al. 2005). 

^ This is equivalent to specifying that the diffusion coefficient 
is proportional to the particle energy. 



and q are the momentum and charge of the accelerating 
particle. Assuming that we are in the strong limit for adi- 
abatic shocks (xc = 4) and that the accelerated charge is 
a highly relativistic electron (/3 = 1, 7 >> 1), we have 

dj 



dt 



1 uqB 2u 

3 TOoC^ c 



By balancing Eq. 4 with Eq. 7 we obtain 



7n 



400uro 




(7) 



(8) 



Eq. 8 is equivalent to Eq. 14 in Eichler & Usov (1993). 
Due to its dependence on _B, 7inax Cb^*- In our standard 
model, 7niax ~ 7 x 10* at the point where the shocks 
intersect the line of centers through the stars. However, 
since the value of C required to match observational data 
is not known a priori, in all of the work presented in this 
paper we fix 7niax = 10^, which is also consistent with the 
maximum 7 at which relativistic electrons are frozen into 
the post-shock fiow. 

2.3. NT emission and absorption 

With a power-law energy distribution, the magneto- 
bremsstrahlung emission can be calculated analytically, 
as in Paper I. However, given an arbitrary energy dis- 
tribution of non-thermal electrons, calculation of the 
magneto-bremsstrahlung emission from the downstream 
flow requires a numerical integration of the magneto- 
bremsstrahlung emissivity over the full energy range of 
the relativistic electron energy spectrum. 



2.3.1. Synchrotron emission 

In the synchrotron limit, the emission arising from a single 
relativistic electron with Lorentz factor 7 is 



P{v) = 5 F{v/vc), 



TOeC^ 



(9) 
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where a is the pitch angle of the particle relative to the 
direction of the B-field, and Fiyjvc) is a dimensionless 
function describing the total power spectrum of the syn- 
chrotron emission, with the frequency where the spec- 
trum cuts off, given by 



37^(7i?sina 3 



47rTOp 



-7 i^bsma. 



(10) 



fb — uj\:,/2t:, where lo\^ = qB/rricC is the cyclotron fre- 
quency of the electron, loq = coh/j is the orbital or rel- 
ativistic gyration frequency of the electron. Values for 
F{u/uc) are tabulated in Ginzburg & Syrovatskii (1965). 
Throughout the remainder of this paper we assume that 
sin a = 1. While this means that our calculations will 
slightly overestimate the synchrotron emission for a given 
C,B, a more detailed consideration is not warranted at the 
present time. At the stagnation point in our "standard 
model", ^ 260 MHz for 7 = 100, and 26 GHz for 
7 = 10^. 

The synchrotron cmissivity formula (Eq. 9) may be 
used even for mildly rclativistic electrons (e.g., 7 < 2), 
provided that x }^ 100, where 



(11) 



(see Fig. 2 of Mahadevan et al. 1996). Since x ^ ^1^-, 
X is small when we have low frequencies and high B- 
fields (or, alternatively when the binary separation is 
small, since B oc l/fgep)- The lowest frequency we exam- 
ine is 200 MHz. In our "standard model" (see Sec. 3.1), 
fb — 1.7 X 10"*^ Hz and B = 5.9 mG at the stagnation 
point of the WCR, which yields x — 1-2 x 10** when 
v = 200 MHz. The declining B-field means that j^b de- 
clines with off-axis distance, which leads to higher val- 
ues of X- the most extreme case which we consider 
(Dsop = 2 X 10^^ cm, Cb = 10"^), x > 100 throughout 
the WCR, justifying our use of the synchrotron cmissivity 
formula. 

2.3.2. Synchrotron self-absorption 

Synchrotron self-absorption (SSA) or, more accurately at 
low 7, magneto-bremsstrahlung self-absorption, is inti- 
mately linked to the energy spectrum of the non-thermal 
particles and their magneto-bremsstrahlung cmissivity, 
and in Appendix A we show how the calculation of SSA is 
performed. The SSA turnover frequency for a source sub- 
tending a solid angle f2 (steradians) is given by (for p = 2; 
cf. Pacholczyk 1970) 



46.3 Bi/65,y^Vi/«17-i/3 Hz, 



(12) 



where B is in Gauss, and Si, is the synchrotron flux in mJy 
at any frequency i> (in Hz) where the source is optically 
thin. For our "standard" GWB model (without IC cooling) 
with C = 10"'^ (see Sec. 3.1), the magnetic field strength 
at the apex of the WCR is 5.9 mG, S^, = 0.026 mJy at 
1/ = 81.1 GHz, and ~ 3 x 10~^^ ster (the spatial extent 



of synchrotron emission from the WCR is ~ 60 x 5 mas), 
giving z/ssA ~ 20 MHz^. For wide CWBs, SSA is negligible 
at GHz frequencies, but may become important in closer 
systems. 

2.3.3. The Razin effect 

When rclativistic charges are surrounded by a plasma 
(as opposed to existing in a vacuum), the beaming effect 
that characterizes synchrotron radiation is suppressed. As 
noted in Paper I, the refractive index of the medium re- 
duces the Lorcntz factor of the electron to 



7 



7 



(13) 



where the plasma frequency i/q = ^JTp^n^Trrrif.. In Paper I, 
this effect was approximated as an exponential reduction 
in flux at the characteristic cut-off frequency. 



v^, = 20 



B 



(14) 



However, since we calculate the magneto-bremsstrahlung 
cmissivity for individual values of 7 in this work, we simply 
replace 7 with 7' in the calculations of the non-thermal 
emission and absorption^. 

2.4. Clumping 

There is a great deal of evidence that the winds of hot stars 
are clumpy (e.g., Moffat et al. 1988; Robert 1994; Lepine 
et al. 2000), and the effect of inhomogeneous winds on 
the WCR is expected to be sensitive to a number of vari- 
ables (e.g., Cherepashchuk 1990; Lepine 1995). Foremost 
amongst these is whether the densities are low enough that 
the WCR and the shocked clumps are adiabatic. This is 
indeed the case in the wide binary systems which wc con- 
sider in this paper*. In such a situation the clumps are 



Note that the calculations in Paper I greatly overestimate 
the strength of SSA (see, e.g.. Fig. 4 in that paper). This error 
resulted from the use of Eq. 6.53 in Rybicki & Lightman (1979) 
that takes C from Eq. 6.20a. We had assumed it was taken from 
Eq. 6.20b, since this is the case for Eq. 6.36. This led to the er- 
roneous conclusion in Paper I that the SSA turnover frequency 
and the characteristic cut-off frequency from the Razin effect 
are comparable in our "standard" CWB model when — 10~^. 



see also the 



We now find that this is the case when C, 
discussion in Sec. 3.5. 

We note that the Razin effect is slightly more subtle than 
replacing 7 with 7' in Eqs. 9 and 10, since not all of the Lorentz 
factors in the function F{v/vc) are related to the beaming 
effect - see, e.g., Rybicki & Lightman (1979), Ginzburg & 
Syrovatskii (1965), or Van Loo et al. (2004) for more details. 
However, this has little impact on the results. 

® The interaction of a clumpy wind with a radiative WCR 
has been discussed by Walder & Folini (2002). In this case 
the WCR may become relatively narrow, though dense, and 
prone to instabilities, such as nonlinear thin-shell instabilities 
(e.g., Stevens et al. 1992; Vishniac 1994). The fate of clumps 
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stripped of material by a combination of dynamical pro- 
cesses (such as ablation) and evaporation (resulting from 
the conduction of heat into their interiors) as they move 
within the WCR. 

The destruction timescale of non-radiative clumps by 
the action of dynamical instabilities is = etcc (Klein 
et al. 1994), where tec is the cloud crushing timescale, 
defined as tec ~ ^c/vs, rc is the radius of the clump, and Vg 
is the shock velocity (for the stationary shocks in CWBs, 
Vs is equal to the pre-shock wind speed), e « 3.5 for a 
density ratio between the clump and the inter-clump wind 
of order 10-100 (Klein et al. 1994), though larger density 
ratios may occur in intrinsic wind shocks. A clump which 
passes through a shock experiences less deceleration than 
the inter-clump material, so a lower limit on the timescale 
for the clump to reach the contact discontinuity is ted = 
Ar/vs, where Ar is the distance from the shock to the 
CD. Clumps will be destroyed before they reach the CD 
when 

id/tcd = rfcc/icd = eJ'c/Ar < 1. (15) 

There are indications that rc ~ li?* at distances of 
~ 20i?* from the star (Lcpine et al. 2000; Rodrigues & 
Magalhaes 2000). Such clumps cannot expand faster than 
their sound speed, unless they behave isothermally, in 
which case they may expand a few times faster. The wind 
temperature of hot stars is ~ 1 — 10 kK, and the sound 
speed of the clump is ~ lOkms^^ . In WR147, the flow 
time of the WR wind to the stagnation point is ^ lO^s, 
so the clumps will have expanded to radii of ~ 1000 R©. 
In our model of WR147 in Sec. 4, Ar « 0.05L'sop, when 
measured along the line of centres of the stars. We thus de- 
termine that td/icd 0.5, which implies that the clumps 
are rapidly destroyed. Material originating in the clumps 
is mixed in with the post-shock gas and the high sound 
speed of the hot WCR works to minimize the density and 
pressure gradients existing in the vicinity of the clump. 
The post-shock flow should rapidly become smooth, and 
the global position of the shocks and the CD should be 
affected very little by the clumpy nature of the winds. 

If td/icd >> 1, clumps could in theory pass com- 
pletely through the WCR and into the unshocked wind 
of the companion star, where they would be surrounded 
by a bowshock and subject to a high level of ablation, 
before being destroyed. The possibility of clumps passing 
through the WCR and interacting with the photosphere 
of the companion star is considered by Marchenko (2002). 

in this situation is less certain. It is possible that they will be 
destroyed when they encounter the WCR, but since the post- 
shock gas is itself cold and dense, new "clumps" may be created 
as the WCR fragments under the action of these instabilities. 
Of note is the fact that the growth rate of the thin-shell in- 
stability depends on high tangential velocities channelling ma- 
terial to the extremities of the kinks (Blondin & Marks 1996) 
- this channelling may be disrupted if the incident flows are 
clumpy, and the growth of this instability moderated to some 
extent, though a numerical calculation would be needed for 
confirmation of this idea. 



Another consideration is whether new clumps are cre- 
ated when gas in a WCR which is largely adiabatic eventu- 
ally cools. We make the general statement that the growth 
of density perturbations within the WCR will depend on 
their size relative to the local cooling length of the sur- 
rounding gas (see Blondin & CiofR 1989), and note that 
any non-thermal component of the pressure (e.g., from 
relativistic particles) lessens the growth rate of the local 
thermal instability by decreasing the dependence of the to- 
tal pressure on the gas pressure (e.g., Wagner et al. 2005). 
Since we know that there are non-thermal particles within 
the WCR we assume that the cold component of the WCR 
is smooth. 

To summarize, we can distinguish the following three 
regions where there may or may not be clumps: stellar 
winds; hot WCR; cold WCR. We have demonstrated that 
clumps are rapidly destroyed in systems where the WCR 
is largely adiabatic, and have noted that clumps are less 
likely to be created within the WCR when there is sup- 
porting pressure from non-thermal particles. Therefore, 
we assume that all parts of the WCR are smooth, leav- 
ing only the nature of the stellar winds to be specified. 
Since the predominant effect of clumping is to introduce a 
scaling factor to the thermal emission, in our "standard" 
model in Sec. 3.1 we assume that the stellar winds are 
smooth. However, our models for WR 147 in Sec. 4 adopt 
clumpy winds in order to match the observed free-free 
emission with specific mass-loss rates. 

2.5. Summary of the calculations 

Our calculation of the observed radio emission from CWB 
systems is based on the density and temperature distribu- 
tion obtained from a hydrodynamical simulation of the 
stellar winds and their collision. Both the stellar winds 
and the WCR contribute thermal emission and absorp- 
tion. Free-free emission and absorption from thermal elec- 
trons is calculated as noted in Paper I, and clumping may 
or may not be included in the model. In addition, the 
WCR has emission and absorption components arising 
from a distribution of non-thermal electrons accelerated 
at the shocks and advected with the downstream flow. 

To describe the non-thermal electron energy distribu- 
tion we use 100, logarithmically spaced, energy bins, ex- 
tending from 7 = 1 to 7max, and which we specify as a 
power-law at the shocks. The pre-cooled normalization of 
this energy spectrum in each cell within the WCR is ob- 
tained by setting the total non-thermal electron energy 
density to a specified fraction of the thermal energy den- 
sity (see Eq. 2). Scalars in the hydrodynamical code are 
then used to apply IC and ionic cooling, and a "cooled" 
non-thermal electron energy distribution is obtained in 
each hydrodynamical cell. As a result, the energy distri- 
bution in each cell within the WCR shows some degree 
of evolution away from its initial power-law (see, e.g.. 
Fig. 3a). In the downstream fiow, the energy loss from 
the non-thermal electrons means that C/rei < C^^th- Since 
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Fig. 2. A density plot of the standard model (showing only the 
positive z-region). The WR star is at (0,0). To aid subsequent 
discussion, two points in the post-shock flow are highlighted: 
a position close to the stagnation point (position "A"), and 
a position some distance ofT-axis in the shocked gas (position 
"B"). 



the magnetic energy density in each hydrodynamical cell 
remains tied to Uth, the magnetic field and non-thermal 
electron energy densities gradually diverge, with Ub ex- 
ceeding Urei in the downstream flow. Electrons which cool 
to 7 = 1 are assumed to return to the thermal pool. 

With the non-thermal electron energy distribution and 
the B-ficld specified in every hydrodynamical cell within 
the WCR, the non-thermal emission and absorption may 
be calculated. If the Razin effect is specified, the de- 
beaming is first calculated and the effective 7 of each en- 
ergy bin is reduced according to Eq. 13. The magncto- 
bremsstrahlung emissivity per electron is then calculated 
using Eq. 9. The magneto-bremsstrahlung emissivity per 
unit volume is obtained from multiplying the emissivity 
per electron by the number density of non-thermal elec- 
trons. Integration over 7 then gives the total non-thermal 
emission per unit volume. Synchrotron self-absorption is 
calculated in a similar fashion. Finally, the radiative trans- 
fer equation is solved to generate synthetic images and 
spectra. 

2. 6. Some simplifying assumptions 

The simplifying assumptions mentioned in Sec. 2.6 of 
Paper I apply also to the work in this paper. An additional 
assumption which was not discussed previously relates to 
the issue of thermal cooling of hot plasma in the WCR. 
This is of relevance since the magnetic energy density Ub , 
and the energy density of relativistic electrons Uici, are 
assumed to be related to the thermal energy density Uth- 
When the separation between the stars is sufficiently small 
that significant cooling of the thermal plasma takes place 
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{E oc 1/Dsep)j a two-fluid implementation of the post- 
shock flow is ideally needed. 

Finally, we emphasize that in our current model the hy- 
drodynamics of the WCR is solved in the limit that a neg- 
ligible amount of energy is placed into non-thermal parti- 
cles. In reality, efficient particle acceleration may greatly 
modify the structure of the shocks bounding the WCR 
and the post-shock flow. Our work in Sec. 4 indicates that 
such changes may actually take place. 

3. Parameter study 

3.1. A standard model 

As with our initial investigations in Paper I, we examine 
the radio emission from a "standard" CWB model with 
the following parameters: Mwr = 2 x 10~^ M0yr~^, 
Mo = 2x10-6 Mgyr-i ,Woo,WR =Woo,o = 2000 kms'^ , 
Z^sop = 2 X 10""^^ cm. The wind momentum ratio, 1] = 
MoVoo.o/ M^!i/^Vao,wR — 0.1, and the distance of the 
stagnation point from the WR and O star is respectively 
''WR = 0.76Dscp and ro = 0.24Z3sop- With such parame- 
ters the WCR is largely adiabatic. We assume solar abun- 
dances for the O star and WC-type abundances for the 
WR star (mass fractions X = 0,Y ^ 0.75, Z = 0.25). 
Temperatures of 10,000 K and an ionization structure of 
H"'", He"*" and CNO^"*" are assumed for the unshockcd stel- 
lar winds. The wind temperatures are lower than those 
used in Paper I, but as the thermal flux from the un- 
shockcd winds is only weakly dependent on the wind tem- 
perature (through the gaunt factor), this change has only a 
small effect on the free- free emission (see Sec. 3.4). As be- 
fore, we adopt a distance of 1.0 kpc for our model system, 
and assume that C = 10"'*. As noted in Sec. 2.3.2, SSA 
is negligible at GHz frequencies in the standard model. 
However, since B (x C^^^i and Si, oc C''/^ in the optically 
thin limit, we find that I'sSA oc The SSA turnover 

frequency is 0.2 GHz when ( = 10^^. 

Since IC cooling is now included in the models, ap- 
propriate luminosities for the stars must also be defined. 
We set iwR = 2 X 10^ and Lq = 5 x 10^ L© for this 
first investigation which results in a photon energy den- 
sity Uph = 2.3 erg cm"-^ at the stagnation point^. For the 
initial non-thermal energy spectrum at the shocks, we as- 
sume a power-law energy distribution with p = 2 and fix 
7min = 1 and 7max = 10^. This latter value is roughly the 
7 at which the energy gain due to shock acceleration is bal- 
anced by the energy loss due to IC cooling, for electrons 
near the apex of the WCR - see Sec. 2.2.3. 

All of our 2D axisymmetric hydrodynamical calcu- 
lations for the standard model use a square grid with 
< r < 4 X 10^5 cm, < z < 4 X 10^^ cm. The z-axis 
is extended to —4 x 10^^ cm with appropriate density and 

® We estimate Lwr and Lq from the following empirical 
formulae. Lwr is calculated from the relationship logAf = 
1.381ogL - 12.1 for Galactic WC stars (Crowther 2003), while 
Lo is determined from logM = 1.691ogL — 15.4 (Howarth & 
Prinja 1989). 
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Fig. 3. The effect of IC and ionic cooling on the relativistic electron energy spectrum (panels a and b) and on the resulting 
synchrotron emission (panels c and d) for a position close to the stagnation point (position "A" , left panels) and some distance 
off-axis in the shocked gas (position "B", right panels), as indicated in Fig. 2. In a) and b) the solid line shows the p = 2 
power-law spectrum assuming no cooling, while the dashed line shows the effects of IC (high energy cut-off) and ionic cooling 
(low energy roll-off). In c) and d) the magneto-bremsstrahlung emissivity resulting from these energy distributions is shown. 
The low frequency cutoff can be influenced by a number of processes. Here, the dashed curves show the Razin effect as an 
example. The emissivity when IC and ionic cooling are or are not included is shown as a solid and dotted curve respectively. 
The difference between the solid and dotted curves indicates the improvements in this paper relative to Paper I. 



temperature values for the unshocked WR wind before 
we process the grid through our radiative transfer code. 
In our models, a viewing angle of 0° implies lines-of-sight 
perpendicular to the axis of symmetry of the model. 



3.2. The thermal flux and r = 1 surface of a hot star 
wind 



To aid the interpretation of the results from our model, we 
calculate the theoretical thermal flux from a single star, 
and the radius of optical depth unity for a sightline from 
infinity. As discussed by Wright & Barlow (1975), the ther- 



mal radio flux for a radially symmetric, isothermal wind 
moving at constant velocity is given by 



4/3 

= 2:^.2 (^] I'JVj'zt^^^mJy, (16) 



where M is in units of MQyr~^ , Vqo is in kms^^ , the 
distance to the system, Z?, is in kpc, and v is in Hz. Zi is 
the rms charge of the ions, 7o is the ratio of electron to 
ion number density (i.e. rie = 7c?^i), Mi is the mean mass 
per ion (in a.m.u.), and ge is the free-free Gaunt factor. 
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Frequency (GHz) 

Fig. 4. Spectra from our standard model with a 0° viewing 
angle - free- free flux (dashed), synchrotron flux (dotted), and 
total flux (solid) . We do not include IC cooling, the Razin effect 
or SSA. With IC cooling included (dot-dash) a clear reduction 
in the amount of synchrotron emission is seen. 

Table 1. Theoretical values of the thermal radio flux, Ss, and 
the radius of the t = 1 surface, Ri, for the winds in our stan- 
dard model. A distance of 1 kpc is assumed, and the winds are 
assumed to be smooth. 



Parameter 


u (GHz) 


WR 





Ss (mjy) 


1.6 


0.322 


0.072 




5.0 


0.637 


0.143 




15.0 


1.216 


0.274 


Ri (10" cm) 


1.6 


4.07 


1.93 




5.0 


1.83 


0.87 




15.0 


0.84 


0.40 



The radius at which a given optical depth from infinity, r, 
occurs, is given by 



Rr = 1.75 X 10 



1/3 1/3^2/3 
28 7 c' V Zj' 

yl/2^1/3 



M 



2/3 



fliVoo V 



cm, 



(17) 



where T is in K. i?o.242 is equivalent to the "characteristic 
radius" defined by Wright & Barlow (1975). 

If the wind is clumped (see also discussion in Sec. 2.4) 
the thermal emission from the star mimics an increased 
mass- loss rate, M = Mactuai/"\/7, where / is the vol- 
mire filling factor of the clumps (Lamers & Waters 1984). 
Values of the rms charge of the ions (^i), the ratio of elec- 
tron to ion number density (70), and the mean mass per 
ion (/ii) for the standard model are 1.12, 1.09 and 4.88 
for the WR star, and 1.00, 1.00 and 1.29 for the O star. 
Values of Se-^ and Rr=i = i?i are given in Table 1. 

3.3. The salient features of the standard model 

Figs. 2 and 3 illustrate the salient features of our mod- 
elling. The density distribution of the unshocked stellar 
winds and WCR in our standard model is shown in Fig. 2. 



Two locations in the post-shock flow are highlighted, one 
near the stagnation point ("A"), and another ("B") at an 
off-axis position. Comparison of Figs. 3a and 3b yields sev- 
eral points worth noting. First, because the non-thermal 
electron energy density at "A" is greater than that at "B" , 
the "initial" energy spectrum (dotted line) has a higher 
normalization in Fig. 3a. Second, the effects of IC and 
ionic cooling are more significant at "A" (solid lines). IC 
cooling is high for flow close to the stars since the inten- 
sity of the radiation field is high, and near the stagna- 
tion point the flow speed is also low, so the exposure time 
is longer. In the standard model, IC cooling is responsi- 
ble for the decline in 7inax from 10^ (the assumed initial 
value) to approximately 200 for flow near the stagnation 
point (Fig. 3a), while for the off-axis position this decline 
stops at ~ 4000 (Fig. 3b). The turn down in the elec- 
tron spectrum at low 7 seen in Figs. 3a and 3b is due to 
ionic cooling. This effect is again strongest for flow near 
the stagnation point due to the higher ionic density and 
slower flow speeds which exist there. 

The magneto-bremsstrahlung emission resulting from 
the "cooled" relativistic energy spectra in Figs. 3a and 3b 
are shown by the solid lines in Figs. 3c and 3d respectively. 
The higher energy electrons that reside at position "B" 
produce magneto-bremsstrahlung emission extending to 
higher frequencies (~ 10^ GHz, versus ~ 10 GHz at posi- 
tion "A" ; see also Eqs. 9 and 10). Thus, flow positions sub- 
ject to high IC cooling produce magneto-bremsstrahlung 
emission over a smaller range in frequency than positions 
with less cooling. With the Razin effect included, a low- 
frequency cut-off to the magneto-bremsstrahlung spectra 
is seen as a result of its de-beaming influence, as shown by 
the dashed curves in Figs. 3c and 3d. It is also worth bear- 
ing in mind that the lower normalization of the relativis- 
tic energy spectrum at "B" (Fig. 3b) versus "A" (Fig. 3a) 
means that at 1 GHz the emissivity per unit volume is 
over an order of magnitude higher at position "A" . 

3.4. Spectral and intensity distributions 

The spectra and intensity distributions obtained from our 
standard model are shown in Fig. 4 and the left column 
of Fig. 5, respectively. For simplicity, the only absorp- 
tion process included in Figs. 4 and 5 is free- free absorp- 
tion. Synchrotron emission dominates the total flux below 
1 GHz. The turnover seen near 400 MHz is due to free- 
free absorption from the unshocked stellar winds. Above 
400 MHz, the synchrotron spectrum is a power-law with 
a spectral index^° a = — (p — l)/2 = —0.5 for p = 2 (see 
Sec. 2). The thermal spectrum is also a power-law with a 
spectral index of -1-0.6, as expected. The slight departure 
of the thermal emission from a power-law at low frequen- 
cies is due to the finite size of our hydrodynamical grid. 
The thermal flux is essentially the sum of the theoretical 



The flux Si, at frequency u, S,y oc u" , where a is the spectral 
index. 
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Fig. 5. The impact of IC cooling on the intensity distributions of our standard model, for a viewing angle of 0° and at 1.6 GHz 
(top) and 22 GHz (bottom). Neither the Razin effect nor SSA are included in this calculation. The images on the left do not 
include IC cooling, while those on the right do. Each image has the same intensity scale and contours. 



flux from each wind (cf. Table 1). However, this is not 
always the case (see Sec. 3.7 and Figs. 10 and 11). 

The effect of IC cooling on spectra from our standard 
model is shown by the dot-dash line in Fig. 4. Its effect 
is greatest at high frequencies, as demonstrated by the 
marked decrease of the synchrotron spectral component. 
However, while this effect is most severe at high frequen- 
cies, IC cooling reduces the flux in the entire 0.1-100 GHz 
frequency range. It is interesting that we do not see a break 
in the synchrotron spectrum, in spite of its clear signature 
in Figs. 3c and 3d. This is because the frequency of the 
break in the synchrotron spectrum is spatially dependent 
(compare Figs. 3c and 3d), and once integrated over the 
entire WCR the frequency break is smeared out and thus 
not visible in the final spectrum. A clearly defined break 
frequency is also lacking in single star models where the 
synchrotron emission arises from shocks intrinsic to the 
wind (Van Loo et al. 2005). 

Spatially, the effect of IC cooling is to reduce the 
brightness of the synchrotron emission from the WCR, but 
more so near the CD than near the shocks (see Fig. 5). 



3.5. Impact of IC cooling and binary separation on the 
synchrotron emission 

The dramatic effect of IC cooling on the resulting spec- 
tra and intensity distributions should be even more pro- 
nounced when the separation between the stars is smaller. 
In Fig. 6 we show the increasing effect of IC cooling on the 
synchrotron emission as the stellar separation is decreased 
from £>sop = 2 X 10^^ cm to Dgcp = 2 x 10^"* cm. Since the 

— 1/2 

intrinsic synchrotron flux scales as Dgcp , we have nor- 
malized the pre-IC cooled spectrum of each model to that 
obtained at Dgcp ~ 2x 10^^ cm. This allows easier compar- 
ison of the divergence between the intrinsic and IC cooled 
spectra as a function of Dsep. The clear trend is for the 
IC cooled synchrotron spectra to steepen with decreasing 

^scp • 

The spectra and intensity distributions from our full 
model (including IC cooling, the Razin effect, SSA and 
free-free absorption) as a function of stellar separation are 
shown in Figs. 7 and 8. While the intrinsic synchrotron 
emission increases with decreasing Dscp, the influence of 
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i^ssA oc C'/'m'^'d-i since cx C^/^M^'ZJ-p^' and 



7/4 1/2 



Sep 



Frequency (GHz) 

Fig. 6. The effect of separation on the IC cooled synchrotron 
spectra with the standard model, a 0° viewing angle and 
( = 10~*, for Daep = 2,5, and 20 x 10^* cm. The spectra are 
normalized so that the intrinsic emission spectra (i.e. without 
free-free absorption, SSA, the Razin effect, or IC cooling) are 
equal to the intrinsic spectrum at -Dscp = 20 x 10^* cm (solid 
line) . 




Frequency (GHz) 

Fig. 7. The effect of separation on the synchrotron spectra 
with the standard model, a 0° viewing angle and = 10~*, 
for Deep = 2, 5, and 20 x 10^* cm. In each case the pre-IC 
cooled intrinsic synchrotron spectrum and the observed IC- 
cooled synchrotron spectrum are shown. 



IC cooling (at high frequencies) and free-free absorption 
and the Razin effect (at low frequencies) means that the 
observed synchrotron emission may decrease with decreas- 
ing separation, which is consistent with the non-detection 
of non-thermal emission from short period CWB systems. 

The low-frequency turnover in Fig. 7 is completely 
due to the Razin effect. Since ric oc AID^^^, and B oc 

(^i/ZjVf^/^^)-! ^ then vji cx (''^/'^M^^^ D'^, and we expect 
the turnover to move to higher frequencies as the separa- 
tion decreases. This is indeed observed, and the turnover 
frequencies arc in excellent agreement with Eq. 14. In con- 
trast, the turnover frequency due to SSA (cf. Eq. 12) scales 



ft oc ZJgcp- Therefore the turnover frequencies resulting 
from the Razin effect and SSA are both oc 1/ Z?sop- If we ig- 
nore free-free absorption for the moment, we find that the 
low-frequency turnover in our standard model is caused by 
the Razin effect rather than SSA when (b ~ Coi ^ 10~^ 
(departures from equipartition are examined in the next 

subsection). Since i^ssa/j^r oc ('^^^M^ ^ , the SSA turnover 
frequency may exceed when M and particularly ^ are 
large. The turnover frequency due to free-free absorption 
depends on the viewing angle into the system, and as it 
has a steeper dependence on Dgep than either the Razin 
effect or SSA (i/ft oc Dso^^'^ - see Paper I), it is increas- 
ingly likely that the low-frequency turnover will be caused 
by free-free absorption as Dgop decreases. 

At the high frequency end of the spectra in Fig. 7, 
the difference between the observed and intrinsic syn- 
chrotron flux is entirely due to IC cooling. With the B- 
field used in the standard model, the emission at a given 
frequency arises from electrons with fairly high values of 7 
(cf. Eq. 10). Thus the influence of IC cooling is fairly high, 
and this is reflected in the curvature in the synchrotron 
spectra at high frequencies. As Dsep decreases, IC cooling 
becomes increasingly important, and we find that high 
frequency synchrotron emission is confined to very thin 
regions immediately behind the shocks. The synchrotron 
emission behind the WR shock can be appreciably higher 
than that behind the O shock, since the non-thermal elec- 
trons are further removed from the O star in this case. 

Fig. 8 shows that the maximum intensity of the emis- 
sion from the WCR may increase with decreasing bi- 
nary separation. However, the bottom panel of Fig. 7 
reveals that the synchrotron emission at 1.6 GHz when 
Dscp = 2 x lO""^* cm is negligible, as it is greatly reduced by 
the Razin effect. The observed emission from the WCR is 
therefore dominated by free-free emission (the increase in 
the free-free emission from the WCR with decreasing sep- 
aration is investigated further in Sec. 3.7). The fact that 
the WCR emission at 1.6 GHz changes from synchrotron 
dominated when Z^scp 2 x 10^^ cm, to free-free domi- 
nated when Dscp = 2 x 10^^ cm, also explains why the "V"- 
shaped emission from the WCR seen at Dgep = 2 x 10^^ cm 
disappears at Dscp = 2 x 10^^ cm (see Fig. 8). 

Fig. 8 also shows that the emission from the WCR 
is sometimes brightest to either side of the apex of the 
WCR, due to high absorption. This is sometimes seen in 
observations (e.g., Watson et al. 2002; Dougherty ct al. 
2005), though it may also arise from time variability in 
the structure of the WCR. 



3.6. Departures from equipartition 

So far, we have assumed that the energy density of the 
non-thermal electrons and the B-field are equal. This as- 
sumption was made for the sake of convenience, but it is 
quite probable that this condition is not attained. For ex- 
ample, the B-field energy density may be equal to the to- 
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mas: 1 



Fig. 8. Intensity distributions at 1.6 GHz and a 0° viewing 
angle with a binary separation of 20 (top), 5 (middle), and 
2 X 10^* cm (bottom), corresponding to angular separations 
of 133, 33 and 13.3 mas respectively. The crosses denote the 
positions of the stars, with the WR star located at (0,0). Each 
image has the same intensity scale and contours. IC cooling, 
the Razin effect, SSA, and free-free absorption are included, 
and = 10~*. Absorption by the WR wind is highest at the 
apex of the WCR - as the binary separation decreases, the 
observed intensity from the WCR reaches maxima either side 
of the symmetry axis. The nature of the emission from the 
WCR also changes with binary separation. In the top panel it 
is predominantly non-thermal, while in the bottom panel it is 
completely thermal. This dramatic change is due to the Razin 
turnover frequency scaling as z/r oc D~^p. 



0.7 

0.6 -'1 

r 

0.5 - 
0.4 - 

I 0.3 - 
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Fig. 9. The synchrotron turnover frequency (normalized to the 
turnover frequency when ("b = (lai) as a function of Cfi/Crei- 
Note that was kept constant in these calculations. We es- 
timate the uncertainty of each data point based on the finite 
frequency sampling in our calculations. The curve is a least 
squares fit of a second order polynomial in logiQ((^B/Crci) to 
logiQ{i/ / vo) , and is forced to go through the position (0.0, 0.0). 
The coefficient's are 0.055, -0.031 and 0.0. 



tal non-thermal energy density from electrons and ions. In 
fact, some astrophysical sources are far from having C/rei = 
Ub (e.g., in some very compact radio sources, Kellermann 
& Pauliny-Toth (1981) note that U.^i/Ub > 10^"!). We 
have therefore explored the effect of varying the ratio 
CB/Crei- In essence, this allows us to explore how changes 
in the B-field affect the emission. The energy density of 
non-thermal ions has no effect on the radio emission. 

When (b is kept constant and d-ei varied, the syn- 
chrotron spectrum simply scales in flux as S*^ cx (^rci (see 
Eq. 6 in Paper I). However, when d-ci is kept constant 
and Cb varied the turnover frequency changes as well as 
the normahzation, since oc 1/B cx At high 

values of Cb we expect SSA to succeed the Razin effect as 
the dominant absorption process behind the low frequency 
turnover in the synchrotron spectrum. When this happens 
we expect the turnover frequency to scale as vssA (x B^^^ 
(see Eq. 12). Our calculations show that it is relatively 
easy to achieve a turnover in the synchrotron spectrum 
at fairly high frequencies through the Razin effect by re- 
ducing Cb, but difficult to achieve turnovers through SSA 
even a little above that obtained when = Coi since the 
dependence of i^ssA on B is weak, and requires commensu- 
rately large values of B. This is illustrated in Fig. 9 where 
we plot the normalized turnover frequency as a function 

of CB/Crol- 

3.7. The thermal flux and binary separation 

The effect of binary separation on the thermal flux from 
the WCR has not previously been considered. As noted 
by Stevens (1995), the free-free cmissivity from the WCR 
is slightly enhanced compared to that from the unshocked 
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Fig. 10. The thermal emission from the standard model with 
a 0° viewing angle and a binary separation of 2 x 10^* cm. 
The thermal emission from the unshocked winds (dashed), 
the WCR (dotted), and the total thermal emission (solid) are 
shown. The intrinsic thermal emission from the WCR (be- 
fore free-free absorption) is shown by the dot-dashed line. 
Synchrotron emission is not displayed. 



winds, at an equivalent distance from the system centre, 
since ef cx p^T^^^^gs. In the WCR the density is typically 
4x higher (cf. the jump conditions for a strong shock), the 
temperature is of order 10^ — 10^ x higher, with perhaps a 
typical increase of a factor of 10"^, and the Gaunt factor is 
approximately double that of the ambient gas. In contrast, 
the free-free absorption is reduced by a factor of ^ lO'^, 
since oc p^T~^^^gg. If the unshocked winds are clumpy, 
and the gas in the WCR basically smooth (see Sec. 2.4), 
the difference in absorption can be even larger. 

In Fig. 10 we show the relative contribution to the 
total thermal emission of the unshocked winds and the 
WCR. Notably, the high frequency thermal spectrum of 
the WCR has a spectral index which is indicative of opti- 
cally thin thermal emission, quite unlike that from the 
unshocked stellar wind. The intrinsic thermal emission 
from the WCR is strongly absorbed at low frequencies, 
mostly by free-free absorption from the unshocked winds 
(i.e. from gas external to the WCR), with a small contri- 
bution from absorption internal to the WCR. The slight 
curvature in the unshocked wind emission at low frequen- 
cies is a result of the finite size of the computational grid. 

We find that the thermal emission from the WCR in- 
creases as the binary separation decreases, while the ther- 
mal emission from the unshocked winds remains broadly 
the same. Since oc cx D~^p and the volume of the 
WCR scales as D^^^, the thermal emission from the WCR 
scales as D~^p. Since the emission from the WCR is op- 
tically thin, we see a substantial fraction of this emis- 
sion at certain viewing angles, and at small separations 
the contribution of the WCR to the total thermal flux 
can be significant. A dramatic consequence of this is that 
a composite-like spectrum (Fig. 10), often presumed to 
be evidence for non-thermal emission, can result entirely 



from thermal processes! The spectral index between 5 and 
15 GHz of the total emission in Fig. 10 is « 0.4, with the 
result that the free-free emission from this model would 
be classified as a "composite" spectrum if it was included 
in the observational analysis of Dougherty & Williams 
(2000). Nevertheless, we believe that a non-thermal origin 
is the most likely explanation in all of the systems noted 
in Dougherty & Williams (2000) since these are typically 
wide systems. 

In Fig. 11a we show how the total thermal flux varies 
with binary separation and viewing angle at = 15 GHz. 
At wide separations and a 0° viewing angle, the total ther- 
mal flux from the system (the sum of the thermal emission 
from the unshocked winds and from the WCR) approaches 
the sum of the theoretical free-free flux from each stellar 
wind. This is because the WCR occurs well outside the i?i 
radius of emission of each wind (see Table 1), and its con- 
tribution to the emission is negligible since the intrinsic 
thermal emission from the WCR scales as D~p. In closer 
systems, the WCR moves within the Ri radius of each 
wind and removes a segment of their emission. However, 
the resulting reduction in flux from the unshocked winds 
is more than compensated for by the increase in ther- 
mal emission from the WCR. Therefore, the total thermal 
emission can considerably exceed that expected from the 
two stellar winds. For example, when Dsep = 2 x lO^^cm, 
the thermal flux from the standard model exceeds the to- 
tal theoretical thermal flux by more than 20%. The ex- 
cess could be considerably higher at smaller separations 
and lower frequencies. The effect is smaller at higher fre- 
quencies, since Ri oc z/""-^. Fig. 11a also shows that the 
observed flux is reduced when one star is in front of the 
other (±90°), as one would expect. 

The nature and degree of clumping may also affect 
how the free-free flux scales with binary separation. If the 
unshocked winds are clumpy, but the WCR smooth due 
to rapid destruction of the clumps, the thermal flux may 
actually decrease with decreasing separation, as shown in 
Fig. lib. This is because as the stars move closer together, 
the increased emission from the WCR is more than offset 
by the decrease in the overall amount of clumpy material 
in the system. Whether the thermal flux increases with de- 
creasing binary separation in a particular system depends 
on the nature and degree of clumping in that system. 

3.8. Synchrotron emission and wind-momentum ratio 

Free-free absorption by the surrounding stellar winds can 
have a large affect on the observable low frequency syn- 
chrotron emission. In previous models where the non- 
thermal emission is treated as arising from a point source, 
the free-free absorption by the stellar winds is determined 
along a single line of sight. Since the spatial extent of the 
WCR is ignored in such models, the low opacity within 
the WCR is also disregarded. In our models, the spatial 
extent of the WCR means that there are multiple lines- 
of-sight to the WCR. The amount of free-free absorption 
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Fig. 11. a) The thermal flux from our standard model at 15 GHz and various viewing angles as a function of the binary 
separation. The dotted line shows the theoretical flux from the WR star, while the dashed line shows the theoretical flux from 
the WR and O star combined. The open circles are for a viewing angle of —90° (WR star in front), the open squares are for 
0°, and the filled squares are for 90° (O star in front). The points at 2 x 10^^ cm and 5 x 10^'* cm are plotted as lower limits, 
since not all of the emission from the WCR is captured in the calculations, b) As a), but for clumpy winds with a filling factor 
of 0.2. In both panels the open circle at 2 x 10^^ cm is hidden underneath the filled square. 



Table 2. Half-opening angles {6, in degrees, measured from 
the line of centres between the stars in our models), of the 
shocks in the WR and O winds bounding the WCR, and of the 
contact discontinuity, as a function of r/. 



V 


WR 


CD 


O 


0.3160 


88 


69 


55 


0.1000 


74 


52 


36 


0.0316 


56 


36 
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that the synchrotron emission suffers varies with viewing 
angle, but is also sensitive to the wind momentum ratio, t], 
since it is the combination of these parameters that deter- 
mines how lines of sight into the system are aligned with 
the opening angle of the WCR. 

In Fig. 12a we show how the synchrotron flux at 
1.6 GHz varies with viewing angle and rj when the bi- 
nary separation is 2 x 10"'^^ cm. The half-opening angle 
of the WR shock, the CD, and the O shock as a function 
of T] are noted in Table 2^^. The main characteristic is an 
increase in absorption of the synchrotron emission when 
the lines of sight to the apex of the WCR pass close to 
either star. 

Here we vary rj by changing Mq, so it is not surpris- 
ing to see that the absorption at -1-90° is highest when rj 
is highest (i.e. Mq highest). When the O star is in front, 
the amount of free-free absorption depends on the ratio of 
R^/ro (Eichlcr & Usov 1993). At 1.6 GHz, this ratio is 

Eichler & Usov (1993) provide an analytical expression for 
the half-opening angle of the CD, 6cd (their Eq. 3). However, 
the statement that the angle between the two shocks, A9 ~ 9, 
is inaccurate, as the data in Table 2 shows. Instead we find that 
A6 ~ 35 — 40° for the range of 77 noted in Table 2, whereas 9 
changes by almost a factor of 2. 



0.58, 0.40, and 0.30 for ri = 0.316, 0.1, and 0.0316, respec- 
tively. Larger values of i?P /ro result in higher absorption, 
which is in agreement with Fig. 12a. We naively expected 
that the synchrotron emission would be least affected by 
absorption when the observer looks down the WCR into 
the system. However, we find that in all 3 cases, the ab- 
sorption is a minimum at viewing angles which are close 
to the asymptotic angle of the WR shock, and steadily 
increases as the viewing angle swings from being parallel 
to the WR shock to being parallel to the O shock. We 
attribute this finding to the fact that even when the view- 
ing angle is within the WCR, the far side of the WCR is 
absorbed by the O star wind, and note that there are im- 
plications for any observing strategy attempting to detect 
radio bursts in CWB systems (e.g., in 7^ Velorum). 

It is also interesting that the model with 77 = 0.316 
shows the lowest relative absorption when the viewing an- 
gle is —90°. Wc attribute this to the wider opening angle 
of the WCR in the model, and the fact that more of the 
WCR is visible when rj is high. Though not shown here, 
the variation with viewing angle is reduced at higher fre- 
quencies, since Ri oc v~^ '^. 

The variation in synchrotron flux with viewing angle is 
more pronounced at smaller binary separations, as shown 
in Fig. 12b, where the values of R^ /tq are 4x larger than 
previously. In all of the models in Fig. 12b, the drop in 
emission as the sightlincs move out of the WCR is steeper 
on the WR side. 

If an observer's line of sight into the system lies in the 
orbital plane, two bursts of radio emission will be seen as 
the WCR sweeps past, as suggested by Eichler & Usov 
(1993). We find that this occurs if R^ /tq > 0.4. When 
R'l/fo ^ 2, the burst is confined to a range in viewing 
angle equal to A9. However, if the inclination of the orbital 
plane is such that the line of sight intercepts the limb of 
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Fig. 12. a) The change in the synchrotron flux at 1.6 GHz as a function of viewing angle using the standard model parameters 
with = 10"'*, but where Mo is adjusted to vary rj. Shown are curves for rj = 0.0316 (solid), rj — 0.1 (dashed), and rj = 0.316 
(dotted). Each curve is normalized to the intrinsic synchrotron emission in each model. SSA, the Razin effect, and IC cooling 
are included. The binary separation is 2 x 10*^ cm. At a viewing angle of —90°, the WR star is in front; at 0° the lines of sight 
are perpendicular to the axis of symmetry of the model; and at +90° the O star is in front. Values for -Rf/ro are noted next to 
each curve. The angular range of the WCR for the model with 77 = 0.316 is shown by the shaded region, with the asymptotic 
angles of the WR and O shocks delimiting its edges, b) As a) but where the binary separation (hence ro) is reduced by a factor 
of 4 to 5 X lO*'* cm. 



the WCR, only one such burst per orbit may be seen. 
In contrast, the behaviour of the system with 77 = 0.0316 
shown in Fig. 12a is better characterized as a single decline 
in synchrotron emission when the WR star is in front of 
the O star, since the emission remains high even when the 
O star is directly in front of the WR star. 

The complexity demonstrated in Fig. 12 further illus- 
trates the need for models which take into account the 
spatial extent of the emission and absorption. 

4. The radio spectrum of WR147 

Having explored how the radio flux is affected by IC cool- 
ing, binary separation, the wind momentum ratio, and 
variations in the B-fleld, we once again turn our attention 
to the modelling of a specific system. WR 147 is one of 
the brightest CWB systems, and is also one of the sys- 
tems where the stellar wind of the WR star and the WCR 
have been resolved at radio wavelengths (Williams et al. 
1997). In addition, WR147 has one of the most exten- 
sively observed radio spectra of any of the CWB systems, 
with observations from 353 MHz to 43 GHz. 

In Paper I, we demonstrated a number of fits of our 
CWB models to the synchrotron spectrum of WR 147. 
However, those models suffer from two main problems. 
Firstly, the optically-thick turnover at ~ 1 GHz was at- 
tributed to SSA which we have discovered subsequently 
was miscalculated (see earlier footnote), and a new fit 
with the corrected model code is now required. Secondly, 
the synchrotron spectrum we were attempting to fit in 
Paper I had been derived from the observed total fluxes 
and an extrapolated thermal spectrum for the stellar wind 
emission based on the flux of the WR star stellar wind at 
several high frequencies (e.g., Setia Gunawan et al. 2001). 



Though a reasonable approach, the resulting synchrotron 
spectrum exhibits considerable variation dependent on the 
adopted thermal model. In an attempt to address this 
problem, here we aim to fit the observed total radio emis- 
sion at each of the observed frequencies, and show the 
contributing synchrotron and thermal spectra from the 
model. 

Observations of WR 147 are available from a number of 
authors. These are summarised in Table 3. Skinner et al. 
(1999) used two methods to obtain flux estimates from the 
same data, and these show considerable scatter. We have 
taken the liberty of quoting an average of the values from 
that paper. It appears that WR147 is variable at some 
frequencies (Setia Gunawan et al. 2001), so where obser- 
vations at multiple epochs are available (e.g., Churchwell 
et al. 1992; Contreras et al. 1996; Contreras & Rodriguez 
1999), we have quoted an average flux. To attain a rea- 
sonable estimate of the uncertainty in the flux values, we 
quote either the variance of the observations (where there 
are multiple observations), or an estimate of the uncer- 
tainty in the absolute flux scale, whichever is largest. For 
the uncertainty in the flux scale, we adopt 5% for fre- 
quencies 15 GHz and greater, and 2% for frequencies less 
than 15 GHz (Perley & Taylor 2003). We believe that our 
approach to estimating the observed fluxes from WR 147 
is very conservative, but it at least makes an attempt at 
assessing the relative merits of all observations in the lit- 
erature, and takes into account the impact of variability. 

WR 147 is a very wide system with a projected sepa- 
ration Dscpcos i = 0.635 ± 0.020"(Williams et al. 1997), 
where i is the angle at which we view the system. At the 
estimated distance of ^ 0.65 kpc (Churchwell et al. 1992; 
Morris et al. 2000) this corresponds to a separation .Dsop ~ 
415/ cos i AU. As noted in Paper I, this relationship 
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Table 3. Observations of WR 147 from the literature. 



Freq. 


Skinner" 


SetiaGunawan'' 


Contreras*^ 


Churchwell" 


Average 


(GHz) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


0.353 




16.0 ±4.0 






16.0 ±4.0 


1.42 


25.2 ± 0.4 


26.4 ±2.0 




26.1 ±0.1^ 


25.9 ± 2.0 


4.86 


37.3 ± 1.6 


35.4 ±0.4 


38.4 ±0.1 


37.5 ±1.0^ 


37.2 ± 1.6 


8.3 


42.5 ± 2.2 




41.0 ±3.7^ 




41.8 ±3.7 


14.94 


47.9 ± 1.7 




57.4 ± 0.3 


47.4 ± 3. 2^ 


50.9 ± 3.2 


22.46 


54.7 ±2.7 






59.1 ± 1.2^ 


56.9 ± 2.8 


43.0 


82.8 ± 1.2 




82.8 ± 1.1 




82.8 ±4.1 



References: a) Skinner et al. (1999), b) Setia Gunawan et al. (2001), c) Contreras et al. (1996), and Contreras & Rodriguez 

(1999) d) Churchwell et al. (1992) 

Notes: 

^ Skinner et al. (1999) quote two values for flux using two different techniques on the same data. The value given here is the 
average value, along with the standard deviation. 

Error based on the standard deviation of multiple observations. 
^ Error assumed to be the uncertainty in the images. 



between -Dgcp and i represents an important constraint 
for any models of the system. We use the stellar wind 
parameters from Paper I: Mwns = 2 x 10~^Moyr~^, 
'I'ooWNs = 950 km s^^ , WqoOB = 1000 km s^^ and Mqb = 
3.8 X lO~^M0yr^^ , giving a wind momentum ratio rj = 
0.02 (Pittard et al. 2002). Of these parameters, A/qb is 
the most uncertain. This is partly due to the consider- 
able uncertainty in the companion's spectral type - es- 
timates range from BO. 5 (Williams et al. 1997) to 05-7 
I-II (Lcpinc et al. 2001). The luminosity of the companion 
(Lob) ranges from ~ 10^ — 10^ Lq, resulting in consid- 
erable scope for varying the degree of IC cooling in the 
models. The composition of the WN8 stellar wind is taken 
from Morris et al. (2000), giving X = 0.09, Y = 0.89, and 
Z = 0.016. The wind temperature for both stellar winds 
was assumed to be 10 kK, and the dominant ionization 
states were H+, He+ and CN02+. 

Firstly, we wish to repeat the modelling done in Paper 
I to evaluate the impact of the miscalculated SSA. The 
resulting models are shown in Fig. 13. As mentioned ear- 
lier, SSA has little impact at the frequencies of the ob- 
servations in these systems, and we have to turn to other 
mechanisms to generate the optically thick turnover. With 
Cb = Crei, we find that the low frequency turnover is 
controlled by free-free absorption along the line-of-sight 
through the circum-binary stellar wind envelope, which is 
contrary to the conclusion from Paper I. The line-of-sight 
opacity changes with the viewing angle of the model: pos- 
itive and negative viewing angles pass first through the O 
star wind and the denser WR star wind respectively. As 
a result, large negative angles are more opaque than large 
positive values, resulting in lower observed fluxes at low 
frequencies with negative viewing angles. The quality of 
the fits in Fig. 13 is largely determined by the model emis- 
sion at 353 MHz, with little variation between models at 
all other observed frequencies. The —30° model matches 
the 353 MHz data the best, though we have made no at- 
tempt to adjust parameters in the model to improve any of 
the fits. Rather, models A-C are intended to demonstrate 
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Fig. 13. Model spectra for WR 147, including total (solid), 
synchrotron (dashed), and thermal (dotted) flux. The obser- 
vations (solid circles) are from Table 3. The open circles are 
synchroton fluxes estimated by subtracting the thermal model 
from the total observed fluxes. Viewing angles of —30° (model 
A), 0° (model B), and 30° (model C) are shown. All models 
used rj = 0.02, 7^ax = 10^ Lqb = 10^ Lq, C = 7.66 x 10-^ 
and a wind filling factor, /, of 0.139. The spectral index of the 
non-thermal electrons, p, was assumed to be 2. A slight cur- 
vature to the optically thin part of the synchrotron spectrum 
(cf. Fig. 4) arises from IC cooling. 

the impact of changing the viewing angle on the behaviour 
of the low frequency spectrum. Cb and Croi are held con- 
stant in models A-C, and their ratio is set to 1. The intrin- 
sic synchrotron luminosity from model B is slightly higher 
than that in models A and C because the stars in model 
B are closer together and this model has a viewing angle 
of 0°. The fit to the 353 MHz data point is very sensitive 
to the viewing angle, and once this point is fit moderately 
well, the goodness-of-fit of the spectral models is largely 
determined by the flux computed at 1.4 and 5 GHz, given 
the smaller uncertainties we have derived for the observed 
fluxes at these two frequencies. We have calculated further 
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models at various viewing angles, and find that we cannot 
obtain satisfactory fits to both the 1.4 and 5 GHz data 
points with = Crci- 

In an attempt to improve the spectral fit, we have cal- 
culated numerous models designed to explore the relevant 
parameter space. Firstly, we varied the wind momentum 
ratio in the range 0.02-0.012 (Pittard ct al. 2002). Since 
the half-opening angle of the WCR is only weakly depen- 
dent on T] in this range {9cd oc rj^^^, see Eichler & Usov 
1993), for a given viewing angle there is very little dif- 
ference to the shape of the spectrum. The normalization 
of the free-free and synchrotron emission changes slightly, 
but this can be compensated for by adjusting the degree 
of clumping and the value of respectively. We also per- 
formed calculations with an increased amount of clumping 
in the undisturbed winds, with the mass- loss rates of both 
stars adjusted downwards in order to obtain the desired 
amount of free-free emission. Changing the volume filling 
factor of the clumps from / = 0.139 to / = 0.03 produces 
an insignificant change in the flux at 353 MHz. 

To simulate crudely a decrease in acceleration effi- 
ciency at the shocks away from the line-of-centers between 
the stars, we have also calculated models where the syn- 
chrotron emission and absorption is limited to a region 
within 7rro/2 of the axis of symmetry of the WCR. While 
this distance is somewhat arbitrary, we note that this 
change has only a small effect on the synchrotron spec- 
trum once it is renormalized (by increasing C). 

The basic problem with all of the models noted so far 
is their total emission spectrum is too flat between 1.5 
and 5 GHz (see Fig 13). The fit in this frequency range 
can be improved by shifting the low-frequency turnover 
in the model synchrotron spectrum to a higher frequency. 
As demonstrated in Fig. 13, fits to these two data points 
cannot be achieved effectively with free- free opacity alone. 
The turnover can be moved to a higher frequency by de- 
creasing the magnetic energy density relative to the rela- 
tivistic energy density, which increases the impact of the 
Razin effect (recall z/r cx 1/B). We induce this reduction 
in the magnetic energy density by setting Croi, as de- 
scribed in Sec. 3.6. The resulting synchrotron spectra are 
shown in Fig. 14, where it is clear that the turnover due 
to the Razin effect moves to higher frequencies as Cb / Crei 
decreases. The best fit to the 1.4 and 5 GHz data is at- 
tained with CB/Crei = 10"'', which also gives the best fit 
to all other data points except that at 353 MHz, which 
is poor. No attempt has been made to improve the fits in 
models D-H by adjusting any other parameters, though 
^B and Crei were adjusted to rcnormalizc the spectrum in 
order to keep the comparison simple. When we reduce Cb 
the synchrotron luminosity drops, and we must increase 
^rei in order to maintain a specific synchrotron luminos- 
ity. This has the consequence that the value of Crci is very 
close to unity in model H. 

A better fit to the 353 MHz data point can be at- 
tained, in part, by increasing the viewing angle to reduce 
the impact of free-free opacity. We have attempted to fit 
the data by varying the viewing angle, Cb and Crci- The 
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Fig. 14. Model synchrotron spectra for WR 147 at a viewing 
angle of 0° and Cb/CicI = 1, and 10"* (models 

D-H respectively, with 10"^ solid - see Table 4). Data points 
are the same as in Fig 13. 

Table 4. Summary of parameters used in our models of 
WR147 with p = 2 (models A-I) and p = 1.4 (models J and 
K). a best fit to the data was only attempted with models 
I, J and K - models A-H are used to demonstrate the ef- 
fect of varying the viewing angle or Cb/CoI. All models in- 
clude free-free absorption, SSA, the Razin effect, and IC cool- 
ing, and were calculated with rj = 0.02 and 7max ~ 10^. 
Models A-I had / = 0.139 and Lob = 10^ Lq, model J had 
/ = 0.170 and Lqb = 10^ Lq, and model K had / = 0.157 
and Lob =4x10^ L© . 



Model Viewing angle (°) Cb/CcI Cd 



A 


-30 


1 


7.66 X 10-^ 


B 





1 


7.66 X 10-^ 


C 


30 


1 


7.66 X 10-=^ 


D 





1 


7.27 X 10-^ 


E 





10-1 


2.00 X 10-2 


F 





10-2 


5.69 X 10-2 


G 





10-^ 


1.72 X 10-1 


H 





10-* 


5.57 X 10^1 


I 


30 


10-* 


5.56 X 10^1 


J 


-30 


1 


1.38 X 10-2 


K 


-60 


1 


1.54 X 10-2 



resulting best fit is shown in Fig. 15 (model I), and its 
weighted value is 10.7. With 4 degrees of freedom, the 
reduced ~ 2.7. At the high frequency end the model 
can be made to more closely fit the data points by reduc- 
ing 7maxj though given the uncertainty in the data points, 
this has little impact on the goodness of the model fit. 

A fundamental concern is the models suggest the DSA 
process in WR147 is highly efficient, with an uncomfort- 
ably large fraction of the pre-shock energy transferred to 
the non-thermal electrons. While there are several ways 
in which the model assumptions might lead us to over- 
estimate C,rci (e.g., the magnetic field in the WCR may 
not be completely tangled resulting in anisotropic emis- 
sion; the thermal energy density of the WCR may have 
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Fig. 15. The total (solid) and synchrotron (dashed) spectra 
obtained with model I (see Table 4), where C^s/Crci ~ 10"'*. 
Also shown is model A for comparison. The thermal emission 
model is the same as used in Fig. 13. 

been underestimated through our choice of mass- loss rates 
and/or wind speeds, etc.), these possible solutions are ei- 
ther rather unsatisfactory in themselves, and/or unlikely 
to substantially increase Cs/Crci- We note instead that in 
situations where the particle acceleration is highly effi- 
cient, the shock precursor (which influences the subshock 
compression, the injection and acceleration eflrciencies, 
and the shape of the non-thermal particle energy spec- 
trum) is modified by the back-pressure of energetic par- 
ticles (e.g., Elhson et al. 2004). As noted in Sec. 2.2, this 
can create a concave curvature to the energy spectrum of 
the non-thermal particles, where p decreases with increas- 
ing energy. Two other effects can change the overall slope 
of the energy spectrum. First, the softer equation of state 
of the relativistic particles makes the shock plasma more 
compressible, and increases the overall compression ratio, 
rtot (when the partial pressure from relativistic particles 
dominates the total pressure of the plasma, rtot — *■ 7). 
Since p = (rtot + 2) /(rtot — 1), a flatter power-law index 
results (p = 1.5 in the limit that rtot = 7). Second, the 
overall compression ratio can become arbitrarily large as 
the highest energy particles escape from the shocks (see 
Ellison et al. 2004, for further discussion), and in this limit 
p^l. 

Since our modelling suggests that the non-thermal par- 
ticle acceleration in WR147 is highly efficient, we have 
explored fits to the data where p is reduced from its "stan- 
dard" value of 2. With <^ioi = Cb and the model parameters 
noted earlier in this section, better fits to the total flux can 
be obtained with p = 1.4. The resulting flt (model J) is 
shown in Fig. 16. C = 1.4 x 10~^ indicates that the particle 
acceleration is indeed efficient, justifying our adjustment 
of p. A p = 1.4 spectrum corresponds to an overall com- 
pression ratio of 8.5 (rtot = (p + 2)/(p — 1)), and is com- 
parable to that achieved in simulations of supernova rem- 
nants with highly efficient particle acceleration (Ellison 
et al. 2004). The high frequency synchrotron fluxes in the 
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Fig. 16. The total (solid), and synchrotron (dashed) spec- 
tra for WR147 with p = 1.4 and 77 = 0.02. Model J has 
Cb = Croi = 1.4 X 10"^ / = 0.17, Lob = lO'' Lq and a view- 
ing angle of —30°. The weighted = 4.97, and the reduced 

= 1.24. Model K has Cb = Croi = 1-54 x 10"^ / = 0.16, 
Lob = 4 X 10^ Lq and a viewing angle of —60°. The weighted 

~ 5.25, but there is only 2 degrees of freedom and the re- 
duced x^ = 2.63. The open circles (triangles) are synchrotron 
fluxes estimated by subtracting the thermal fluxes in model J 
(model K) from the total observed fluxes. The synchrotron 
data points are different between models J and K and from 
those in Figs. 13-14 because of the changes in the thermal flux 
by way of the filling factor. 



model can be reduced by increasing Lqb, while maintain- 
ing a similar quahty of flt (model K, see Fig. 16), though 
we again emphasize that the flts to the WR 147 data are 
to the total flux, not to model deduced synchrotron fluxes 
which are very uncertain. 

An important point is that as the overall compres- 
sion increases (rtot > 4), the subshock compression ra- 
tio declines. Since the subshock is responsible for heat- 
ing the thermal gas, the post-shock temperature drops 
below the value expected in the strong-shock ideal-gas 
limit. Analysis of the X-ray data from WR147 indicates 
that the temperature of peak emission from the thermal 
gas is kT « 0.9 keV (Skinner et al. 1999). This can be 
compared to the theoretical value, kT = 1.96/iw|keV, 
expected after thermalization in the strong-shock ideal- 
gas limit for material of mean particle mass fimn flowing 
into a stationary shock at a typical wind terminal veloc- 
ity Woo = lOOOuskms"^ . As most of the X-ray emission 
should arise from the shocked WN8 wind (see Pittard & 
Stevens 2002), we set /i = 1.17 and = 0.95. The ex- 
pected temperature is then kT = 2.1 keV. While the con- 
trast is not as large as that seen in WR 140, where the 
theoretical temperature is 5 times higher than what is ob- 
served (Pollock et al. 2005), we nevertheless conclude that 
there is perhaps some additional support for highly effi- 
cient particle acceleration. 
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5. Summary and future directions 

We have extended the work in Paper I by including the 
evolution of the non-thermal electron energy distribution 
in the post-shock flow. IC cooling is by far the most im- 
portant process to consider, but we have also calculated 
the degree of ionic cooling for completeness. 

We have demonstrated that the geometry of the WCR 
and the streamlines of the flow within it lead to a spa- 
tially dependent break frequency in the synchrotron emis- 
sion. We therefore do not observe a single, sharp break in 
the synchrotron spectrum integrated over the WCR, but 
rather a steepening of the synchrotron spectrum towards 
higher frequencies. 

As expected, IC cooling reduces the overall brightness 
of the synchrotron emission from the WCR, but more so 
near the CD than near the shocks. The synchrotron emis- 
sion appears brighter at the WR shock than at the O shock 
when the O star dominates the luminosity of the system. 
In this situation, IC cooling has a somewhat lesser effect at 
the WR shock. Current radio interferometers do not have 
the spatial resolution and sensitivity needed to detect such 
spatial variation, but the predicted morphology may allow 
future observations to determine whether the non-thermal 
electrons are primarily accelerated at the shocks through 
the DSA scheme. 

We find that the observed synchrotron emission may 
decrease with decreasing Dgcp, through a combination of 
IC cooling, free-free absorption, and the Razin effect. This 
is consistent with the non-detection of non-thermal emis- 
sion from short period CWB systems. 

We have also investigated how the synchrotron emis- 
sion varies with the assumed B-field. and have discovered 
that it is very difficult for SSA to produce high turnover 
frequencies (> 1 GHz). It is much more likely that the 
cause of these turnovers is free-free absorption through 
the circumbinary envelope, or the Razin effect. If the op- 
tically thin thermal radio emission from the WCR, which 
scales as -D~p, becomes comparable to the thermal emis- 
sion from the unshocked winds, it is possible to obtain a 
composite-like spectrum entirely from thermal processes. 

Furthermore, the total free-free emission from CWB 
systems is likely to vary with binary separation, though 
the exact nature of the variation depends on the degree 
and distribution of clumping in the system. This has im- 
portant consequences for deriving mass-loss rates for stars 
in binary systems using radio fluxes, with the implication 
that they may be significantly over- or under-estimated. 
Our models show that it is preferable to use a high fre- 
quency (> 15 GHz) for estimates of mass-loss rates. 

We have applied our new model to WR 147, and find 
that the introduction of IC cooling allows a better fit 
to the high frequency data than obtained in Paper I. 
However, since a weighted measure is biased by the rel- 
atively low uncertainity of the 1.4 and 5 GHz data points, 
any model which fits these points will be preferred, and 
in this respect, the other data points are almost irrelevant 
given their relative weights. With the data currently avail- 



able, we find that the best fits to the 1.4 and 5 GHz fluxes 
are obtained when either the magnetic field energy density 
is substantially less than the non-thermal electron energy 
density, or when the injected non-thermal electron spec- 
trum has a power-law index of p = 1.4. We find that we 
require efficient electron acceleration in both cases, which 
leads us to favour the model where p = 1.4, since non- 
linear calculations of DSA show a hardening of the non- 
thermal particle spectrum. 

In summary, our analysis of the radio data from 
WR147 indicates that efficient particle acceleration oc- 
curs in this system, and that the shock structure may 
be modified by the back-pressure of energetic particles. 
We emphasize that our current models are in many ways 
very simplistic, and in order that these conclusions can be 
given a firmer footing, more sophisticated models which 
self-consistently include the back-reaction of the acceler- 
ated particles on the shock structure need to be developed. 
Such calculations, unfortunately, are not trivial, and a fur- 
ther complication is that the shock obliquity may change 
with off-axis distance. As the obliquity increases, the in- 
jection rates may reduce, so that some parts of the shocks 
may be highly modified, while other parts will be unmod- 
ified. 

Such self-consistent models would also be of more gen- 
eral use, because despite the apparent efficiency of DSA, 
X-ray spectra from CWBs have been modelled and in- 
terpreted assuming that the shocks place an insignificant 
fraction of their energy into non-thermal particles (e.g., 
Stevens ct al. 1996; Zhekov & Skinner 2000; Pittard & 
Corcoran 2002). These works have led to inferences for 
important quantities such as the mass-loss rates, wind 
speeds and composition, and the rate of electron and pro- 
ton equilibration. However, the strong coupling between 
particle acceleration and thermal heating implies that the 
inferences made from X-ray observations may differ sub- 
stantially between interpretations that include particle ac- 
celeration self-consistently and those that do not. Future 
work should address this issue and models should attempt 
to fit radio and X-ray data simultaneously. Finally, to fur- 
ther improve our understanding of particle acceleration in 
WR147, and other colliding wind binaries, accurate and 
precise radiometry is required in order to distinguish be- 
tween models. 
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Appendix A: Synchrotron self-absorption (SSA) 

For an isotropic electron distribution function f{p), the 
absorption coefficient for SSA is given by (see Rybicki & 
Lightman 1979, Eq. 6.46): 

a. = J 47r/ dp ifipi) - f{p)) P., (A.l) 
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where P,^ is the magneto-bremsstrahlung emissivity and 
pI = (7^ - l){mecf + [hv/cf - 2hi^-fme. (A.2) 



Using/(p) = 0f , wcfind that 



n{Ei) n{E) 
S-Khiy^ J \ piEi pE 

where 

El ^ E — hv = jnipC^ — hv. 



P^pE dE, 



(A.3) 



(A.4) 



Re-writing the integral in Eq. A.3 in terms of 7, we get 




piEi pE 



Pvl\/l'^ - 1 d7. (A.5) 



We evaluate Eq. A.5 to find the SSA for any isotropic 
distribution of particles 71(7). 
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